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ABSTRACT 

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave 
signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to 
measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 
2015. We study an astrophysically motivated population of sources (binary components with masses 1.2 Mq- 
1.6 Mq and spins of less than 0.05) using the full LIGO analysis pipeline. While this simulated population 
covers the observed range of potential binary neutron-star sources, we do not exclude the possibility of sources 
with parameters outside these ranges; given the existing uncertainty in distributions of mass and spin, it is critical 
that analyses account for the full range of possible mass and spin conhgurations. We find that conservative prior 
assumptions on neutron-star mass and spin lead to average fractional uncertainties in component masses of 
^ 16%, with little constraint on spins (the median 90% upper limit on the spin of the more massive component is 
^ 0.7). Stronger prior constraints on neutron-star spins can further constrain mass estimates, but only marginally. 
However, we find that the sky position and luminosity distance for these sources are not influenced by the 
inclusion of spin; therefore, if LIGO detects a low-spin population of BNS sources, less computationally 
expensive results calculated neglecting spin will be sufficient for guiding electromagnetic follow-up. 

Subject headings: gravitational waves - methods: data analysis - stars: neutron - surveys 


1. INTRODUCTION 

As we enter the advanced-detector era of ground-based 
gravitational-wave (GW) astronomy, it is critical that we un¬ 
derstand the abilities and limitations of the analyses we are 
prepared to conduct. Of the many predicted sources of GWs, 
binary neutron star (BNS) coalescences are paramount; their 
progenitors have been directly observed |Lorimer| ( |2008 1 , and 
the advanced detectors will be sensitive to th eir GW emission 
up to ^ 400 Mpc away (jAbbott et al.|2016a i. 

When analyzing a Gw signal frorn a circularized compact 
binary merger, strong degeneracies exist between parame¬ 
ters describing the binary (e.g., distance and inclination). To 
properly estimate any particular parameter(s) of interest, the 
marginal distribution is estimated by integrating the joint pos¬ 
terior probability density function (PDF) over all other parame¬ 
ters. In this work we sample the posterior PDF u sing software 
implemented in the LALInference library (|Veitch et al. 


farr@uchicago.edu 

^ Enrico Fermi Institute and Kavli Institute for Cosmological Physics, 
University of Chicago, Chicago, IL 60637, USA 

^ School of Physics & Astronomy, University of Birmingham, Birming¬ 
ham, B15 2TT, UK 

^ Canadian Institute for Theoretical Astrophysics, 60 St. George Street, 
University of Toronto, Toronto, Ontario, M5S 3H8, Canada 

^ Department of Physics, University of Maryland-College Park, College 
Park, MD 20742, USA 

^ Gravitational Astrophysics Lab, NASA Goddard Space Flight Center, 
Greenbelt, MD 20771, USA 

^ The Pennsylvania State University, University Pai'k, PA 16802, USA 

^ Leonard E. Parker Center for Gravitation, Cosmology, and Astro¬ 
physics, University of Wisconsin-Milwaukee, Milwaukee, WI 53201, USA 

^ LIGO Laboratory, California Institute of Technology, Pasadena, CA 
91125, USA 

^ Astrophysics Science Division, NASA Goddard Space Flight Center, 
Code 661, Greenbelt, MD 20771, USA 

Massachusetts Institute of Technology, 185 Albany St, Cambridge, 
MA 02139, USA 


2015|), specihcally we use results from LALInfernc e_nest 
( Veitch and Vecchio|2010 1 , a nest samp ling algorithm ( Skilling| 
2UU6|l, and LALlNFERENCE_MCMC ([Christensen et al.|2004[ 



gnds; 


a Markov-chain 
chapter 12). 
largely assessed pa¬ 
rameter constraints assuming negligible neutron-star (NS) spin, 
restricting models to nine parameters. This simplihcation has 
largely been due to computational constraints, but the slow 
spin of NSs in short-period BNS syst ems observed to date 
(e.g.,[Mandd and O’Shaughnessy|2010|l has also been used for 
justification. However, proper ch^acterization of compact bi¬ 
nary sources must account for the possibility of non -negligible 
spin, otherwise parameter esti mates will be biased ( |Buonanno| 
|et al.|2009) [Berry et al.|20f5l l. This bias can potentially lead 
to incorrect conclusions about source properties, and even 
misidentihcation of source classes. 

Numerous studies have looked at the BNS parameter esti¬ 
mation abilities of ground-based GW detectors such as the 
Advanced La ser Interferometer Gravitational-Wave Observa- 
tory (aLIGO;|Aasi et al.|[2015| l and Advanced Virgo (AdV ; 
[Acernese et al.|2015| l detectors. [Nissanke et al.| ( |2()10| 2011| ) 
assessed localiz ation abilities on a simulated non-spinning 
BNS population. |Veitch et al.|(|2012]l looked at several poten¬ 
tial advanced-detector networks and quantlhed the parameter- 
estimation abilities of each ne twork for a signal from a fiducial 
BNS with non-spinning NSs. |Aasi et aL] ( |2013[ l demonstrated 
the ability to characterize signals from non-spinning BNS 
sources with waveform models for spinning sources using 
Bayesian stochastic samplers in the LALI NFERENCE library 
(|Veitch et al. 12015)1. |Hannam et al.|(|20131l used approximate 
methods to quantify the degeneracy between spin and mass es¬ 
timates, assuming the compact objects’ spins are align ed with 
the orbital angular momentum of the binary (but see |Haster| 
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et al.||2015|l. |Rodriguez et al.|(|2014|l simulated a collection 
of loud signals from non-spinning BNS sources in several 
mass bins and quantified parameter estimation capabi lities in 
the advanced-detector era using non-spinning models. |Chatzi-| 
|ioannou et al.|(|2015|l introduced precession from spin-orbit 
coupling and found that the additional richness encoded in 
the waveform could reduce the mass-spin degeneracy, helping 
BNSs to be distinguish ed from NS-black hole (BH) binaries; 
Littenberg et al. (2015|l conducted a similar analysis of a large 


catalog of sources and found that it is difficul t to infer the pres 
ence of a mass gap between NSs and BH s (Ozel et ^2010 


|Farr et ari|2011[ |Kreidberg et al.|[2012 1, although this may 


(Mandel et al. 

2015 

. Finally, Singer et al.| 

follow-on [Berry et a 

.(2015 

1 represent an (al 


and the 


during the first 1-2 years of the advanced-detector era. These 
studies simulated GWs from an astrophysically motivated BNS 
population, then detected and characterized sources using the 
search a nd follow-up tools that are used for L IGO-Virgo data 
analysis (|Aasi et al.|2014[| Abbott et al.|2016c|l. The final stage 
of the analysis missing from these studies is the computation¬ 
ally expensive characterization of sources while accounting 
for the compact objects’ spins and their degeneracies with 
other parameters. The present work is the final step of BNS 
characterization for the|Singer et al.|(|2014)l simulations using 
waveforms that account for the effects of NS spin. 

We begin with a brief introduction to the source catalog used 
for this study and Singer et al. (20141 in section]^ Then, in 
section]^ we describe the results of parameter estim ation from 
a full analysis that includes spin. In section [XT] we look at 
mass e stima tes in more detail, and spin-magnitude estimates in 
section l3T2l In section|4]we conside r the estimation of extrinsic 
parameters: sky position (section [4~T|l and distance (section 
|4.2^ , which we do not expect to be significantly affected by the 
inclusion of spin in the analysis templates. We summarize our 
findings in section A comparison of computational costs 
for spinning and non-spinning parameter estimation is given 
in appendix [A| 

2. SOURCE SIMULATION AND SELECTION 

We have restricted our study to the first year of t he advanced- 
detector era, using the same 250 simulations that|Singer et al. 
(|2014|l analysed with non-spinning parameter estimation. For 
these, Gaussian noise was generated using the ‘early’ 2015 
aLIGO noise curve found in Barsotti and Fritsch^ ( |2012l l. Ap¬ 
proximately 50,000 BNS sour ces were simulated, using t he 
SpinTaylorT4 waveform model ( |Buonanno et al.|2003||2009| l, a 
post-Newtonian inspiral model that includes the effects of pre- 
cession, to generate the GW signals. Component masses were 
uniformly distributed between 1.2 M© an d 1.6 Mf?), which 
reflects the range of observed BNS masses ( Ozel et al.|2012 1 . 
Component spins were isotropically oriented, with magnitudes 
Xi .2 = c\S\^ 2 \/Gm\ 2 drawn uniformly between 0 and 0.05; 
here |S'i. 2 | are the NSs’ spin angular momenta and mi 2 their 
mass, the indices 1 and 2 correspond to the more and less 
massive components of the binary, respectively. The range of 
simulated spin magnitudes was chosen to be consistent with 
the observed population of short-period BNS systems, cur- 
rently bounded by PSR J0737-3039A ( |Burgay et al1[2003 


Brown et al. |2012|l from above. Finally, sources were dis 


duce a netw ork signal-to-noise ratio (S/N) of pnet = 5 (Singer 
et al. 2014 1 , where pnet = X); P? individual detector S/Ns 
Pi combined in quadrature. 

Of this simulated population, detectable sources were se¬ 
lected usi ng the GSTLAL_INSP IRAL matched-filter detection 
pipeline (|Cannon et al. 112012)1 with a single-detector S/N 
threshold p > 4 and false al^m rate (FAR) threshold of 
FAR < 10“^ yr“*. The FAR for real detector noise is largely 
governed by non-stationary noise transients in the data that can 
mimic GWs from compact binary mergers, which Berry et al. 
( |2015| l demonstrate make negligible difference to parameter 
estimation for the (low-FAR, BNS) signals considered here. 
Because our simulated noise is purely stationary and Gaussian 
with no such artifacts, FAR estimates are overly optimistic. 
To compensate, an additional threshold on the network S/N 
of Pnet >12 was applied. This S/N threshold is consistent 
with the above FAR th reshold when applied to data similar to 
previous science runs ( [Abbott et aH2016al|Berry et al.|2015] l. 
A random subsample of 250 detections were selected for pa- 
rameter estimation with LALInference[]3] This mass and 
spin distributions of this subset is sta tistically consisten t with 
those the sources w ere drawn from ( [Berry et al.|2015 i. See 
[Singer et al.[ ( [2014| ) for more details regarding the simulated 
data and GSTLAL_inspiral analyses. 

3. SPINNING ANALYSIS 


[Singer et al.[ ( [2014| l details the detection, low-latency local¬ 
ization, and medium-latency (i.e. non-spinning) follow-up of 
the simulated signals in 2015. In this work we perform the ex¬ 
pensive task of full parameter estima tion that accounts fo r non¬ 
zero compact-object spin. Whereas [Singer et al.[([2014[l used 
the (non-spinning) TaylorE2 wavefor m model, we make use 
of the SpinTaylorT4 waveform model ( [Buonanno et al.|2003[ 
[2009[ l, parameterized by the fifteen parameters that uniquely 
define a circularized compact binary inspiral p] 

We assume the objects to be point masses with no tidal 
interactions. The estimation of tidal parameters using post- 
Newtonian approximations is rife with systematic uncertainties 
that are comparable in magnitude to sta tistical uncertainties 
^Yagi and Yunes[2014[[Wade et al.|2014[ l. Though marginaliz¬ 
ing over uncertainties in tidal parameters can affect estimates 
of other parameters, the fact that tidal interactions only impact 
the evolution of the binary at late times (only having a mea¬ 
surab le impact at frequencies above ~ 450 Hz; Hinderer et al. 
[2010[ ) limits both their measurability and the resulti ng biases 
in other parameter estimates caused by ignoring them[Damour| 

[etaLlISfet . 

The simulated population of BNS systems contains slowly 
spinning NSs, with masses between 1.2 Mq and 1.6 Mq and 
spin magnitudes x < 0.05. This choice was motivated by the 
characteristics of NSs found thus far in Galactic BNS systems 
expected to merge within a Hubble time through GW emission. 
However, NSs outside of B NS systems have been observed 
with spins as high as x = 0.4 ([Hessels et al.|2006[ [Brown et al. 
2012^, and depending on the NS equation of state (EOS) could 


theoretically have spins as high as x ^ 0.7 (Lo and Lin|201 Ijl 


tributed uniformly in volume (i.e. uniform in distance cubed) 
to a maximum distance at which the loudest signal would pro¬ 


** The mean (median) pnet of the set of 250 events is 16.7 (14.6). 

The fifteen parameters are two masses (either component masses or 
the chirp mass and mass ratio); six spin parameters describing the two spins 
(magnitudes and orientations); two coordinates for slcy position; distance; an 
inclination angl e; a polarization ang le; a reference time, and the orbital phase 
at this time (see|Veitch et al.|2015| for more details). The masses and spins 
are intrinsic parameters which control the evolution of the binary, the others 
are extrinsic parameters which describe its orientation and position. 
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without breaking up. For these reasons, the prior assumptions 
used for Bayesian inference of source parameters are more 
broad than the spin range of the simulated source population. 

To simulate a real analysis scenario where the class of com¬ 
pact binary and the NS EOS are not known, we use uniform 
priors in component masses between 0.6 Mq and 5.0 Mq to 
avoid any prior constraints on mass posteriors, and our standard 
BH spin prior: uniform in spin magnitudes X 1.2 ~ 6^(0,1) and 
isotropic in spin orientation. Prior distributions for the location 
and orientation of the binary match that of the simulated popu¬ 
lation, i.e. isotropically oriented and uniform in volume (out to 
a maximum distance of 218.9 Mpc, safely outside the detection 
horizon, which is ~ 137 Mpc for a 1.6 M 0 -I .6 Mq binary) 
Choosing any particular upper bound for spin magnitude woula 
require either assuming hard constraints on NS spin-up, which 
are based upon observations with hard-to-quantify selection 
effects, or making assumptions regarding the unknown EOS 
of NSs. For these reasons we choose not to rule out compact 
objects with high spin a priori by using an upper limit of X< 
encompassing all allowed NS and BH spins. In section [T3] 
we look at more constraining spin priors, and particularly how 
such choices can affect mass estimates. 

We describe parameter-estimation accuracy using several 
different quantities, depending upon the parameter of interest. 


Simplest is the fractional uncertainty {x), where ax 
and {x) are the standard deviation and mean of the distri¬ 
butions for parameter x respectively. This is particularly 
useful for showing how uncertainty scales with S/N: in 
the limit of high S/N, the standard deviation can be ap¬ 
proximated from the (i nverse) Fisher m atrix and scales 
inversely with the S/N ( Vallisneri|2008 1 . 


• The credible interval CI^ in the range that contains 
the central p of the integrated posteri or, with (l-p )/2 
falling both above and below the limits ( |Aasi et al.|2013[ ). 
Specifying the credible interval for several values of p 
gives information about the shape of the posterior. 


• As an alternative to credible intervals, we use credible 
upper or lower bounds. These are the one-sided equiv¬ 
alents of credible intervals, and are useful for distribu¬ 
tions that are peaked towards one end of the parameter 
range or for parameters we are interested in putting a 
limit upon (the spin magnitude satisfies both of these 
criteria). 


• For sky-localization, we use credible regions (the two- 
dimensional generalization of the credible interval) 
which are the smallest sky area that encompasses a given 
total posterior probability. The credible region for a total 
posterior probability p is defined as 

CRp = argmax [ dS7Pa(fl), ( 1 ) 

A Ja 

where Pa,{il) is the posterior PDF o ver sky position H , 
and A is the sky area integrated over ( |Sidery et al.|2014[ ). 
We also consider the searched area A*, the area of the 
smallest credible region that includes the true location. 

To check that differences between our spinning and non¬ 
spinning analyses were a consequence of the inclusion of spin 


The mean (median) true distance for the set of 250 events is 52.1 Mpc 
(47.8 Mpc), and the maximum is 124.8 Mpc. 


and not because of a difference between waveform approxi- 
mants, we also ran SpinTaylorT4 analyses with spins fixed to 
X\. 2 = 0. There were no significant differences in parameter 
estimation between the non-spinning TaylorF2 and zero-spin 
SpinTaylorT4 results for any of the quantities we examined ^ 
Therefore, we only use the TaylorF2 results to illustrate the 
effects of neglecting spin. 


3.1. Mass Estimates 

To maximize sampling efficiency, model parameterizations 
are chosen to minimize degeneracies between parameters. To 
leading order, the post-Newtonian expansion of the wave¬ 
form’s phase evolution depends on the chirp mass. Me = 
(mim 2 )^^^(mi -l-m 2 )“*^^, making it a very well constrained pa¬ 
rameterization of binary mass. The second mass parameter 
used is the mass ratio q = m 2 /mi, where 0 < q < 1. Detectors 
are much less sensitive to the mass ratio, and strong degen- 
eracies with spin make constraints on q even worse (|Cutler| 
|and Flanagan|1994| . It is primarily the uncertainty in q that 
governs the uncertainty in component masses mi and 1112 - 

Figure[2shows the superimposed, one-dimensional marginal 
posterior PDFs and cumulative density functions (CDFs) for 
the chirp mass (centered on each mean) and mass ratio for 
all 250 events. As a representation of a typical event’s poste¬ 
rior distribution, we show the average PDFs and CDFs, where 
the average is taken over all 250 posterior PDFs and CDFs at 
each point. Chirp-mass distributions are usually well approx¬ 
imated by normal distributions about the mean, while mass 
ratio estimates have broad support across most of the prior 
range, whereas the simulated population had a narrower range 
between 0.75 and 1. 

To trace individual parameter uncertainties across the pop¬ 
ulation we use the fractional uncertainties in chirp mass 
aMc! (-^c) and mass ratio Uqj(q). The chirp mass and mass 
ratio conveniently cover mass space (which is why they are 
used for sampling), but the total mass M = m \+1712 is also of in¬ 
terest for determining the end product of the merger, so we also 
plot the fractional uncertainty um! (M). The mean (median) 
fractional uncertainties in chirp mass, mass ratio and total mass 
for the simulated population are 0.0676% (0.0642%), 28.7% 
(28.4%) and 6.15% (5.81%) respectively. For comparison, the 
mean (median) fractional uncertainties in chirp mass, mass ra¬ 
tio and total mass from the non-spinning analysis are 0.0185% 
(0.0165%), 8.90% (8.79%) and 0.542% (0.491%) respectively. 
We furthe r exa mine the impact of spin on mass measurements 
in section 

The fractional uncertainties for the chirp mass, mass ratio 
and total mass all decrease as S/N increases, as shown in Fig¬ 
ure]^ which also shows results from the non-s pinning analysis. 
As expecte d from Fisher-matrix studies (e.g., Finn and Cherj 
|noff||199^ , most appear to be inversely proportional to the 
S/N: the exception is a^j {q) from the spinning analysis, which 


Using as an example the chirp mass, the most precisely inferred parame¬ 
ter, we can compare the effects of switch from a non-spinning to a spinning 
analysis to those from switching waveform approximants by comparing the 
difference the posterior means (Ale). The difference between means from the 
SpinTaylorT4 analyses with and without spin, is an order of magnitude greater 
than the difference between the zero-spin SpinTaylorT4 and TaylorF2 anal¬ 
yses: defining the log-ratio (, = logj^CKAlc}® - Alc}°|/| (Ale)”® - Alc)°|), 
where the superscripts S, 0 and NS indicates results of the fully spinning 
SpinTaylorT4, the zero-spin SpinTaylorT4 and the non-spinning TaylorF2 
analyses respectively, the mean (median) value of ^ is 0.90 (1.04), and 92.4% 
of events have ^ > 0 (indicating that the shift in the mean from introducing 
spin is larger than the shift from switching approximants). 
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better fit as cx Pne/ . We do not suspect there is anything fun- 
damental about the c>c , rather it is a useful rule-of-thumb. 
The behaviour can still be understood from a Fisher-matrix 
perspective, which predicts a Gaussian probability distribution 
(with width oc Pnet)- Since the mass ratio is constrained to be 
0 < ^ < 1, if the width of a Gaussian is large, it is indistinguish¬ 
able from a uniform distribution and the standard deviation 
tends to a constant 1 j'/Vl ~ 0.289. When the width of the 
Gaussian is small (< 0.1), the truncation of the distribution is 
negligible and the standard deviation behaves as expected, as is 
the case for the non-spinning results. The standard deviations 
obtained for the spinning runs lie in the intermediate regime, 
betwe en being independent of S/N and scaling inversely with 
it (cf. ILittenberg et al. 112015b: the mean (median) standard 
deviation Oq is 0.182 (0.183) ^ The mass-spin degeneracy 
broadens the posteriors for botnuie chirp mass, the mass ratio 
and the total mass; a consequence of the broadening for the 
mass ratio is that the uncertainty does not decrease as rapidly 
with S/N (over the range considered here). 

Projecting the tightly constrained chirp mass and poorly con¬ 
strained mass ratio 90% credible region from space into 
component-mass space makes it obvious how important mass- 
ratio uncertainties are for extracting astrophysical information. 
The credible regions in component-mass space are narrow ba¬ 
nanas that lie along lines of constant chirp mass, bounded by 
the constraints on mass ratio (see Figurej^for some examples 
posteriors). 


3.2. Spin Estimates 

We now look at the constraints placed on the spin of the 
slowly spinning simulated BNS sources. Even though the sim¬ 
ulations occupy a small fraction of the spin-magnimde prior 
volume, most posterior distributions span the majority of the 
prior range. For non-precessing systems, where the orbital 
plane is stationary with respect to the line-of-sight, varying the 
spin of the compact objects has a similar effect on the phase 
evolution of the GW as varying the mass ratio. This results in 
a strong degeneracy between the two parameters. Modulation 
of the GWs from precession of the orbital plane can break this 
degen eracy (jVecchio|2004[[Lang and Hughes|2006[[Vitale et ak] 

|2014||Chatziioannou et al.|2015|i; however, only systems with 

large spins that are misaligned with the orbital angular momen¬ 
tum signihcantly precess. Non-precessing systems, with either 
low or aligned spins, suffer the most from this degeneracy as 
the only information regarding the mass and spin is encoded 
in the phase of the GW. The simulated sources in this study 
fall in the latter category of low spins. Figure shows the 
distribution of Gaussian kernel density estimates of the PDFs 
for the spin of the most and least massive components, xi and 
X 2 , respectively. The labeled regions of figure Inbound the 
specihed percent of PDFs as a function of spin, where the 90% 
region, for example, is bounded by the 5th and 95th percentiles 
of the PDFs at each spin value. 

The spin of the more massive component has a larger effect 
on the GW, and is therefore systematically better constrained, 
as seen in Figure For both spins, however, the posterior 
shows slow spins to be only slightly more probable than high 
spins for most sources. The mean (median) 50% upper limits 

The uncertainty for the symmetric mass ratio rj = +m 2 )^, 

which is constrained to be 0 < r; < 1/4, does scale approximately as p/j,. The 
mean (median) standard deviation for the spinning runs is 2.00 X 10“^ 
(1.95 X 10-2). 


on xi and xi are 0.319 (0.302) and 0.424 (0.419) respectively; 
the 90% upper limits are 0.707 (0.699) and 0.855 (0.859). 


3.3. Prior Constraints on Spin 

Since spin is largely degenerate with mass ratio, and spin 
is expected to be small for BNS sources, it is interesting to 
ask how the mass constraints are affected by making stronger 
prior assumptions about the spin of NSs. First, we make the 
extrem e assumption that N Ss have negligible spi n, as was 
done in jSinger et al.| ( |2014| l and jBerry ^ al.|p015| l. Figure]^ 
compares the distribution of (fractional) uncertainties in chirp- 
mass, mass-ratio and total-mass estimates for the spinning and 
non-spinning analyses. The average fractional uncertainties 
from the non-spinning analysis are a factor of ~ 3^ smaller 
than the uncertainties from a spinning analysis in the case of 
Aic and q, and about an order of magnitude smaller for the 
total mass. 

Figure compares cartoon 90% credible re gions in 
component-mass space of 5 chosen simulated signals (|Hannam| 
et al.|201^|Chatziioannou et al.|2015 cf.). As a consequence 
of the difficulty of estimating the narrow and nonlinearly cor¬ 
related credible regions in m\-m 2 space, we illustrate the cred¬ 
ible regions in mi-m 2 space as the projection of a rectangular 
region in Ai^-q space. To define the rectangular region we 
use 90% credible intervals of the one-dimensional posterior 
PDFs of Aic and for Aic we use the central 90% credi¬ 
ble interval (5th to 95th percentile), and for q the upper 90% 
credible interval (10th to 100th percentile). These differing 
credible intervals were chosen to better summarize the one¬ 
dimensional posterior PDFs, which are typically normal for 
Aic and skewed toward high values for q (see hgure[T]). 

We can investigate the impact of stronger prior assumptions 
regarding the maximum spin of NSs on mass estimates by dis¬ 
carding posterior samples above a given spin. Figure|^shows 
the cumulative distribution of lower 90% bounds on me esti¬ 
mates of m 2 among the 250 simulated sources for spin priors 
of Xi, 2 < {1,0.7,0.4,0}. Xi, 2 < 1 and xi ,2 = 0 correspond 
to the spinning and non-spinning analyses described above. 
X < 0.7 is consistent with the NSs remaining intact for most 
proposed non-exotic EOSs. x < 0.4 is consistent with the spin 
of observed, isolated NSs to date. 

Erom these PDEs, it is clear that fairly strong prior assump¬ 
tions on NS spin are required to significantly impact mass 
constraints. Assuming NSs to be spinning with xi. 2 < 0.4 a 
priori only constrains masses by an extra few percent compared 
to allowing them to have Xi, 2 < 1- 


4. SOURCE LOCATION 

Having discussed how GW observations can measure the 
intrinsic properties of their source systems, we now consider 
the measurement of e xtrinsic parameters, specih cally the sky 
position (section |TT] l and the distance (section |4~2) l. These 
are central to the success of multimessenger astronomy. The 
sky position is required in order to direct telescopes for elec¬ 
tromagnetic (EM) follow-up and to verify that any observed 
transients do coincide with the source of the gravitational 
waves. The distance also aids electromagnetic follow-up as it 
allows cross-reference wit h galaxy catalogs to find the most 
probable source locations (jNissanke et al.|2013HHanna et al.j 

|2014[|Fmi et al. 120141 [Blackburn et al.|2015| l. Even wimout an 

observed counterpart, me posterior for me (three-dimensional) 
position allows us to assign a probability that the source re¬ 
sides in given galaxies; combining the redshift of these galaxies 
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Figure 1. Superimposed posterior probability density (top) and cumulative density (bottom) functions for the chirp mass and mass ratio of all spinning analyses. 
The solid lines show the average distribution for the simulated population. The chirp-mass distributions have been centered on the distributions’ means to highlight 
their consistent morphology. 


(measured electromagnetically) with the gravitational-wave 
luminosity distance gives a measure of the Hubble constant 
free of the usual systematics (|Schutz|1986||Del Pozzo|2012|l. 
For our population of slowly spinning NSs, we do not expect 
the measurement of the extrinsic parameters to be affected by 
the inclusion of spin in the analysis. 


4.1. Sky Localization 

In order for EM observatories to follow-up a GW detection, 
they need an accurate sky location. This must be provided 
promptly, while there is still a visible transient. Parameter esti¬ 
mation while accounting for spin is computationally expensive 
and slow to complete (see appendix [A)l. There are alternative 
methods that can provide sky localization more quickly. The 
most expedient is BAYESTAR, which uses output from the de- 
tection pipeline to rapidly compute sky position dSinger and 

ice||2016| l. BAYESTAR can compute sky positions with a 
latency of a few seconds. Between the low-latency BAYESTAR 
and the high-latency full parameter estimation, there is the 
medium-latency option of performing non-spinning parameter 
estimation with computationally cheap TaylorF2 waveforms. 
This requires hours of wall time to complete analyses, with 
the exact time depending upon the degree of parallelization. 
Despite only using information from the detection triggers, 
rather than full waveforms, it has been shown that BAYESTAR 
produces sky areas for BNS signals fully consistent with non¬ 
spinning parameter estimation results, prov ided that there was 
a trigger from all detectors in the network (Singer et al.|2014 


[Berry et al.|2015j [Singer and Price|2016| l! Having now per¬ 

formed a full spinning analysis, we can compare the results 
of high-latency parameter estimation with the more expedient 
methods of inferring sky position. 

In Figure 1^ we show the cumulative distributions of re¬ 
covered 50% credible regions, 90% credible regions, and 
searched areas. All three quantities show good agreement 


across all parameter-estimation techniques. For the for the 
slowly-spinning BNSs considered here, including spin in the 
analysis does not change the average ability to localize sources 
on the sky. 

We can consider sky localization in greater detail by com¬ 
paring areas on an event-by-event basis and not just the cu¬ 
mulative distribution across the population. Doing this, we 
confirm that sky localization is consistent between approaches 
for any given event. We use the medium-latency, non-spinning 
TaylorF2 analysis as a reference point, and compare the ratio of 
sky areas. To summarize the variation in sky areas computed 
in different analyses, we use the log ratio 

- l°gl0 ^^NS^ ’ 

where is a credible region or the searched area as deter¬ 
mined by method X and is the same quantity from the 
non-spinning analysis. The log ratio TZ^ is zero when analysis 
X agrees with the non-spinning results. Considering all 250 
events, the mean and standard deviation of the log ratio is given 
in Tabled For the purposes of EM follow-up, there is no sig¬ 
nificant mfference between analysesp^The computationally 
expensive fully spinning analysis does not improve sky local¬ 
ization: there is no disadvantage in using the lower-latency 
results for EM follow-up of slowly spinning BNSs. 


4.2. Luminosity distance 
The distance is degenerate with the inclination ([Cutler and 


Flanagan|1994[ Aasi et al.|2013 1, and the inclination can be 


better constrained for precessing systems (van der Sluys et al. 


The non-spinning analysis was performed with LALlNFERENCE_NEST 
while the spinning analysis was performed with LALlNFERENCE_MCMC 
(Veitch et al.|20151; therefore the consistency between analyses additionally 
shows the consistency of results from different sampling algorithms. 
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Signal-to-noise ratio 



Figure 2. Fractional uncertainties in chirp mass A4c, mass ratio q and total 
mass M estimates as a function of network S/N for both the fully spinning 
SpinTaylorT4 analysis and the medium-latency non-spinning TaylorF2 analy¬ 
sis. We only show statistical uncertainties, not systematic errors (which are 
present when spin in not included). The lines indicate approximate power-law 

trends (oc for spinning cr^/(q) and oc for the rest) to guide the eye. 


2008 11 Vitale et al. |2014|l . Since we are considering a population 
with low spins, precession is minimal, and there should be little 
effect from including spin in the analysis. 


The absolute size of the distance credible interval CI^ ap¬ 
proximately scales with the distance, hence we divide the 
credible interval by the true (injected) distance D*; this gives 
an approximate an alogue of twice the fractional uncertainty 
( |Berry et al.|2015| l. The cumulative distribution of the scaled 
credible intervals is plotted in Figure The mean (median) 



Figure 3. The distribution of one-dimensional marginalized posterior prob¬ 
ability density functions (PDFs) of spin magnitudes of the more and less 
massive components (xi and xi^ respectively) for all 250 simulated sources. 
Shaded regions show the 90% credible boundaries for the spin distributions 
of the population, and the solid lines show the average of each PDF. The pos¬ 
teriors have consistent morphology and span the majority of the prior range. 
The spin of the most massive component is typically slightly more constrained 
toward low values, but even a maximal spin of xi = 1 is never ruled out with 
100% certainty. 


Table 1 

Comparison of sky localization areas produced by the 
low-latency BAYESTAR analysis and the high-latency fully 
spinning SpinTaylorT4 analysis to those produced by the 
medium-latency non-spinning TaylorF2 analysis. The 
mean and standard deviation of the log ratio for the 50% 
credible region CRq.s, the 90% credible region CR 0.9 and 
the searched area A* are listed for each analysis. 



BAYESTAR 

Spinning 

Logarithmic 

Standard 

Standard 

ratio 

Mean deviation 

Mean deviation 


^^CRo.5 

0.095 

0.117 

0.022 

0.067 


0.075 

0.094 

0.028 

0.063 


0.106 

0.447 

0.002 

0.397 


values of CIq^/D^, for the spinning and non-spinning analy¬ 
ses are 0.436 (0.376) and 0.426 (0.363) respectively; the val¬ 
ues of CI^g/D* are 0.981 (0.845) and 0.951 (0.819), and the 
fractional uncertainties aoj{D) are 0.302 (0.262) and 0.245 
(0.239). There is negligible difference between the spinning 
and non-spinning analyses as expected. 

5. CONCLUSIONS 

In this study we investigated the effects of accounting for 
spin when estimating the parameters of BNS sources with 
aLIGO. We expect NSs to be only slowly spinning, and hence 
that their spins only have a small effect of the GW signature 
of a BNS merger. However, allowing for spins does have a 
significant effect on parameter constraints. Strong degenera¬ 
cies are present in the model; not only are the spins themselves 
poorly constrained, but these degeneracies result in weaker 
constraints on other parameters, particularly masses. Exclud¬ 
ing spin from parameter estimation results in artificially narrow 
and potentially inaccurate posterior distributions. 

Weaker constraints are the result of accounting for broad 
prior assumptions on NS spins. We tested various choices 
for conservative prior assumptions about NS spins and found 
them to have little effect on mass estimates. Only strong 
prior assumptions, such as say x\, 2 0.05 (consistent with 
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Figure 4. Fractional uncertainty in chirp mass j\4c, mass ratio q and total 
mass M estimates from the non-spinning and spinning analyses. The mean 
fractional uncertainties from the non-spinning analysis are 0.0185%, 8.93% 
and 0.542% for chirp mass, mass ratio and total mass respectively. These are a 
factor of a few smaller than found from a spinning analysis (0.0676%, 28.7% 
and 6.15% for chirp mass, mass ratio and total mass respectively). 


the simulated population, and NSs observed in short-period 
BNS binaries to date) are likely to significantly affect mass 
constraints. However, such strict prior assumptions are hard 
to justify given the small number of observed systems and 
possible selection effects. 

We performed parameter estimation on a astrophysically 
motivated population of BNS signals, assuming an aLIGO 
sensitivity comparable to that expected throughout its first 
observing run. Using a prior on spin magnitudes that is uni¬ 
form from 0 to 1, spanning the range permitted for BHs and 



m.2/M0 


Figure 5. Approximate 90% credible regions for the component-mass esti¬ 
mates of 5 selected simulations from the spinning analysis; each region is the 
projection of a rectangular region of chirp-mass-mass-ratio space, bounded 
by the central 90% credible interval in chirp-mass and upper 90% credible 
interval in mass-ratio. Circles indicate the true masses of each simulation, and 
bars indicate the lower bounds of the upper 90% credible intervals (i.e. the 
10th percentiles) on mass ratio for increasingly strict prior assumptions on the 
maximum spin of NSs. 



Figure 6. Cumulative distributions of the lower bounds of the upper 90% 
credible interval (i.e. 10th percentiles) on the estimated mass of the least 
massive binary components under increasingly strict prior assumptions about 
maximum NS spin. Restricting spins to be below break-up (x $ 0.7) for 
non-exotic equations of state has little effect, as does restricting spin to the 
maximum observed NS spin (x $ 0.4). Only strict prior assumptions on NS 
spin will significantly impact mass constraints. 


extending beyond the expected (but uncertain) upper limit for 
NSs, the median 90% upper limit on the spin of the more 
massive component is 0.70 and the limit for the less massive 
component is 0.86. The median fractional uncertainty for the 
mass ratio Uqj{q) is ^ 30%, the median fractional uncertainty 
for the total mass um / (M) is ~ 6% and the median fractional 
uncertainty for the chirp mass crAtc/(-^c) is ~ 0.06%. Despite 
the mass-spin degeneracy and only weak constraints on the 
spin magnitudes, we hnd that we can place precise constraints 
on the chirp mass for these BNS signals. 

The sky-location accuracy, which is central to performing 
EM follow-up, is not affected by including spin in the analysis 
of low-spin systems; this may not be the case when spin is 
higher, i.e. in binaries containing a BH. Eor our population 
of BNSs, sky localization is unchanged by the inclusion (or 
exclusion) of spin in parameter estimation. The median CRo.g 
(CRo.s) is ~ 500 deg^ (^ 130 deg^). The luminosity distance is 
similarly unaffected for this population of slowly spinning NSs; 
the median fractional uncertainty uoj(p) is ~ 25%. However, 
an analysis that includes spins requires the use of more com- 
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Figure 7. Cumulative fractions of events with sky-localization areas smaller 
than the abscissa value. Top: Sky area of 50% credible region CRq.s Middle: 
Sky area of CRo.g. Bottom: Searched area A,. The high-latency results 
including spin are indicated by the solid (blue) line . The lower latency non¬ 
spinning and BAYESTAR from|Singer et al.H2014J are denoted by thicker 
(green) and thinner (orange) lines respectively. The 68% confidence intervals 
for the cumulative distribution are denoted by the shaded areas. 


putationally expensive waveforms (that include more physics), 
increasing latency by an order of magnitude. Therefore, if the 
population matches our current expectation of being slowly 
spinning, the low-latency results that could be supplied in time 
for EM observatories to search for a counter-part are as good as 
the high-latency results in this respect, and there is no benefit 
in waiting. 

Following t he submission of thi s article, aLIGO made its 
first detection ([Abbott et al.|2016bll. This was of a binary BH 


system ([Abbott et al.|2016c|) rather than a BNS, but much of 
our understanding of the abilities of the parameter-estimation 
analysis, such as the effects of mass-spin degeneracy, trans¬ 
lates between sources. The era of GW astronomy has begun, 
and parameter estimation will play a central role in the science 
to come. 
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APPENDIX 


A. COMPUTATIONAL COST 

Performing a fully spinning analysis is computationally ex¬ 
pensive. The main computational cost is generating the Spin- 
TaylorT4 waveform, which must be done each time the like¬ 
lihood is evaluated at a different point in parameter space. 
Progress is being made in reducing the cos t of generating 
waveforms and evaluating the likelihood (e.g.,[Canizares et al. 
2013 [ [Ihirrer|2014J . Employing reduced order modelling can 
speed up the non-spinning TaylorF2 analysis by a factor of 
~ 30 ( [Canizares et al.[20T5] l. This is still to be done for a wave¬ 
form that includes the effects of two unaligned spins; however, 
progress has also been made in constmcting frequency domain 
approximants using shifted uniform asymptotics, which can 
speed up generatio n of a waveform li ke SpinTaylorT4 by an 
order of magnitude ( [Klein et al.|2014[ l. 

In figure we present the approximate wall time taken for 
analyses comparable to those presented here. The low-latency 
BAYESTAR and the high-latency fully spinning SpinTaylorT4 
results are for the 250 events considered here. The medium- 


latency non-spinning TaylorF2 results are from Berry et al. 
(2015 these are not for a different set of signals, but represent 
a similar population (in more realistic non-Gaussian noise). 
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Scaled credible interval Scaled credible interval Clj^g/D* 

Figure 8. Cumulative fractions of events with luminosity-distance credible intervals (divided by the true distance) smaller than the abscissa value. Left: Scaled 
50% credible interval CI^j/D*. Right: Scaled 90% interval CI°g/D*. Results using the spinning analysis are indicated by the solid (blue) line and the results 
using the non-spinning analysis JBerry et al.|2015) are indicated by the dashed (green) line. The 68% confidence intervals for the cumulative distribution are 
denoted by the shaded areas. 


representing what we hope to achieve in reality[^ The wall 
times f or BAYESTAR are s ignificantly rednced compared to 
those in|Berry et al.|(|2015|l because of recent changes to how 
BAYESTAR integrates over distance ( Singer and Price|2016| l: 
the mean (median) time is 4.6 s (4.5 s) and the maxirnum is 
6.6 s. We assume that 2000 (independent) posterior samples 
are collected for both of the LALInference analyses. The 
number of samples determines how well we can characterize 
the posterior: ~ 2 000 is typically needed to calculate CRq.q 
to 10% accuracy ( |Del Pozzo et al. |2016 l. In practice, we 
may want to collect additional samples to ensure our results 
are accurate, but preliminary results could also be released 
when the medium-latency analysis has collected 1000 samples, 
which would after half the time shown here with a maximum 
wall time of 5.87 x 10"^ s ~ 16 hr. We see that the fully spinning 
analysis is significantly (here a factor of ~ 20) more expensive 
than the non-spinning analysis, taking a mean (median) time 
of 1.48 X 10® s ~ 17 days (9.19 x 10® s ~ 11 days) and a 
maximum of 1.48 x 10® s ~ 172 days. 

The times shown in figure [^illustrate the hierarchy of times 
associated with different analyses. However, they should not 
be used as exact benchmarks for times expected during the 
first observing run of aLIGO because the version of LAL¬ 
INFERENCE used here are not the most up-to-date versions. 
Lollowing the detection pipeline identifying a candidate BNS 
signal, we expect BAYESTAR results with latency of a few sec¬ 
onds, non-spinning LALINFERENCE results with a latency of 
a few hours, and fully spinning LALINFERENCE results only 
after weeks of computation. 

While work is underway to improve the latency of and to 
optimize parameter estimation with LALINFERENCE, there is 
also the possibility of developing new alg orithms that provide 
parameter estimates with lower latency (|Haster et aLpOlS" 


Pankow et al.|2015i. Improving computational efficiency is 


important for later observing runs with the advanced-detector 
network: as sensitivities improve and lower frequencies can 
be measured, we need to calculate longer waveforms (at even 
greater expense). 
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